Intermittent dry granular flow in a vertical pipe 
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The intermittent compact flow of glass beads in a vertical glass pipe of small diameter is studied 
experimentally by combining particle fraction, pressure, and air and grain flow rates measurements 
with a spatio-temporal analysis of the flow. At the onset of the flow, a decompaction front is observed 
to propagate from the bottom to the top of the tube at a velocity much larger than that of the 
grains. The blockage front also propagates upwards and at a still higher velocity. The decompaction 
induces a decreasing pressure wave strongly amplified as it propagates upwards towards the top of 
the tube. Pressure variations of 3000 Pa or more are detected in this region while particle fraction 
variations are of the order of 0.02. Grain velocities during the flow period also increase strongly 
at the top of the tube while the corresponding fraction of total time decreases. A ID numerical 
model based on a simple relation between the effective acceleration of the grains and the particle 
fraction variations reproduces the ampliflcation effect and provides predictions for its dependence 
on the permeability of the packing. 

PACS numbers; 45.70.-n, 81.05.Rm 



I. INTRODUCTION 

Dry particle flows driven by gravity (as in emptying 
hoppers or silos) or by a gas current (in pneumatic trans- 
port are encountered in many industrial processes. 
Modeling such flows represents a challenging fundamen- 
tal problem since they involve complex interactions of 
moving grains with the surrounding air, bounding walls 
and other grains. Various flow regimes may be observed 
in these systems according to the spatial distribution of 
the particle concentration and velocity fleld. Examples 
of such regimes include the free-fall of particles at high 
velocities and low particle fractions, slow compact flows 
with high particle fractions, and density waves in which 
compact and dilute zones alternate 0- These particle 
flows may be stationary (in this case, particles - or den- 
sity waves - propagate at a constant mean velocity) or 
display oscillations, stick-slip motions, intermittency or 
blockage effects. Such non-stationarities have been re- 
ported in the simultaneous flow of air and grains in ver- 
tical pipes H, H, IE 111 and are encountered in a variety of 
settings. They represent a signiflcant problem in many 
large-scale industrial facilities such as silos or pneumatic 
transport systems: intermittent flows may induce very 
large, potentially destructive, pressure variations. 

Several approaches have been suggested to model theo- 
retically granular flows in pipes, each of them being appli- 
cable to only a part of the flow regimes. Kinetic theories 
0, 13 are best suited when the particle fraction is lower 
than 0.5 and coUisional effects play an important part. 
In this regime, the time between two successive collisions 
is large compared to the duration of the collision. At 
higher particle fractions, grains are in contact with their 
neighbors a large fraction of the time and friction forces 
between particles and walls become significant. Although 
attempts have been made to adapt the kinetic theory to 
these frictional regimes jOj] , different approaches have of- 
ten been used. In their study of compact moving bed 
flows, Chen et al. ^3] use the classical Janssen theory 



[ill Il2| to estimate friction forces on the walls from the 
relation between the vertical component of the stress- 
tensor in the grain packing and the horizontal normal 
stress component on the walls. 

The present work deals with intermittent compact 
flows of grains in a vertical tube of small diameter. In 
this case, the particle fraction is very close to that of a 
static random packing of particles, and the solid friction 
between particles and walls is always important. The 
objective of this study is to analyze the spatio-temporal 
characteristics of the onset and blockage of such inter- 
mittent flows. Of special interest is the relation between 
air pressure, particle fraction and grain velocities and the 
manner in which variations of these quantities are ampli- 
fied as they propagate along the flow tube. 

A similar flow regime is encountered in the "ticking 
hourglass" experiment 0, 0, 0, 0, 0| correspond- 
ing to a periodic intermittent granular flow between two 
closed glass containers connected by a short vertical con- 
striction. In this case, intermittency is believed to be 
largely due to the build-up of a pressure difference in- 
duced by the grain flow between the two containers (both 
because air is dragged downwards and because the empty 
volume available for air increases with time in the upper 
container and decreases in the lower one). Once the back 
pressure is sufficiently large to halt flow, it decays due 
to air flow through the grain packing in the constriction 
until it is low enough to permit flow. In this experiment, 
the key phenomena are localized in a small volume close 
to the constriction; in particular, static arches of particles 
are believed to build up right above it, where transient 
fractures of the packing are also observed. 

In the intermittent compact flow down a long vertical 
tube considered in the present study, the flow intermit- 
tency phenomenon is distributed over most of the height 
of the tube. This allows one to perform simultaneously a 
large variety of measurements (local pressures and par- 
ticle fractions, air and grain flow rates, . . . ) at several 
locations, and so characterized the variations of the flow 
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FIG. 1: View of the flow experiment. - Insert: sketch of the 
constriction at the bottom end of the pipe. 



with height. High resolution spatio-temporal diagrams 
of particle motions along the tube are also obtained and 
allow for the characterization of propagation of the onset 
and blockage of the grain flow along the tube length. 

In section^ we describe the experimental set-up and 
procedure. Then, in section IIIII the qualitative charac- 
teristics of the compact intermittent flows are described, 
and quantitative data and estimates of air pressure, par- 
ticle fraction and flow rate variations are presented. In 
section Hvl we discuss appropriate volume and mass con- 
servation equations and their bearing on the experimen- 
tal results. The various forces acting on the grains and 
the manner in which they may induce intermittency are 
discussed in section^ The equations of motion are then 
written and solved numerically in order to model the spa- 
tial variations of pressure, particle fraction and grain ve- 
locities. 



II. EXPERIMENTAL SET-UP AND 
PROCEDURE 

The granular flow takes place in a vertical glass tube 
of radius J2= (1.5 ±0.025) mm and length L= 1.25 m 
(Fig. Grains flow from a spherical hopper welded 
to the top of the pipe. The bottom end is fitted with 
a variable constriction allowing one to adjust the out- 
flow: this constriction is realized by a cylinder with 
a conical tip which can be precisely moved across the 
pipe (see insert of Fig. QJ. The grains are spherical 
glass beads of diameter 2a = (175 ± 25) /J,m and density 
p = (2.50 ± 0.02) 10^ kgm^^. Electronic scales are placed 
below the tube and their output is transmitted to a com- 
puter at the rate of one sample per second. The mass 
flow rate Fm of the grains is then determined from the 
time variation of the measured mass. The grain flow rate 



will be characterized by the mean superficial velocity q 
representing the volume flow rate of grains per unit area 
with: 



F„ 



(1) 



The hopper is connected to the tube but otherwise closed 
except for an air inlet. The volume flow rate Q"^'^'^ of air 
into the hopper is measured by an on-line sensor (the 
precision is ±1% at the largest flow rates). An addi- 
tional transducer measures the air pressure inside the 
hopper. Continuity requires that the inflow of air com- 
pensates for the volume of both air and grains leaving 
through the vertical tube. Therefore if ph remains con- 
stant with time, Q'^°'^ is the sum of the volume flow 
rates of air and of the grains in the experimental tube. 
The superficial velocity Qa of air in the tube (volume flow 
rate per unit area) thus satisfies: 
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Four pressure sensors are placed along the tube at vertical 
distances 200, 450, 700 and 950 mm below the outlet of 
the hopper. They are connected to the inside of the tube 
through 0.5 mm diameter holes in the tube wall. A fine- 
meshed grid is stretched across the holes on the external 
wall to prevent grains from entering the sensors. An in- 
dependent experiment has been realized to estimate the 
pressure drop across the grid and the portholes during 
typical pressure transients encountered in the intermit- 
tent flows. Due to the low dead volume of the transducer 
(23 /il), the corresponding error on the pressure measure- 
ment is less than 10 Pa for the fastest transients: this will 
be seen below to be much smaller than the pressure to 
be measured. 

Local mean particle fraction variations in tube sec- 
tions are measured by an electrical capacitance sensor 
using two shielded 3 mm diameter cylindrical electrodes 
pressed against the outside tube walls and facing each 
other. The sensor is connected to a Gi?1620 capacitance 
bridge (10 kHz measurement frequency) and a lock- in de- 
tector with a 680 fis time constant. The noise level is 
of order 10~^ pF, which corresponds to particle fraction 
variations of order 10"'^. 

In order to test the device, an air-grain interface was 
moved through the sensor and the output was recorded 
as a function of the location of the interface. This allowed 
one to determine that the measurement is averaged over 
a slice of the sample of typical thickness 2 mm parallel to 
the tube axis. The variation of the DC output with the 
particle fraction is assumed to be linear and the probes 
are calibrated by comparing readings obtained with an 
empty tube and with the same tube filled with a static 
bead packing. This assumption of a linear variation is 
frequently used in the literature: in the present case, it is 
particularly justified by the small particle fraction varia- 
tions during the experiment and by the random, approx- 
imately isotropic structure of the packing. The particle 
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fraction Cmax of the static bead packing obtained after 
filling the tube is assumed to be equal to 0.63: this value 
corresponds to a Random Close Packing of spheres and 
is frequently obtained experimentally in similar systems. 
A direct measurement of Cmax has been realized but its 
precision is limited by the small amount of beads in the 
tube. A slightly lower value Cmaa; = 0.61 ± 0.02 was ob- 
tained (this difference may be due to the small diameter 
of the tube increasing the influence of wall effects reduc- 
ing locally the particle fraction). In view of the large 
uncertainty on this measurement, the value 0.63 is re- 
tained in the following. 

The capacitance measurements can be performed at 
chosen heights by moving the sensors along the tube. 
The instantaneous outputs of four pressure and capaci- 
tance sensors are digitized simultaneously at a sampling 
frequency of up to 1 kHz by a Spectral Dynamics 51)195 
signal analyzer. Time-averaged values of all pressure, ca- 
pacitance and flow rates sensors are also recorded but at 
1 s intervals. 

The dynamical properties of the flow are studied 
from spatio-temporal diagrams constructed from appro- 
priately assembling the output of a digital linear CCD 
camera. Light intensity variations along a vertical line 
precisely aligned with the tube are recorded at sampling 
rates of 500 lines per second with a resolution of up to 
2048 pixels (the equivalent pixel size on the tube ranges 
from 10 /im to 500 /xm depending on the magnification). 
One thus obtains 2D images in which the x-axis repre- 
sents time while the z-axis corresponds to distance along 
the tube. 



III. EXPERIMENTAL RESULTS ON THE 
INTERMITTENT COMPACT FLOW 

The flow regimes are controlled by adjusting the con- 
striction at the bottom of the experimental tube. The 
compact flow regime is observed for narrow constrictions, 
at typical superficial velocities q between 0.01 ms~^ and 
0.08 ms~^. Narrowing further the bottom constriction 
leads to a blockage of the flow while density waves ap- 
pear at higher flow rates. The typical aperture 6 of the 
constriction (see Fig. ^ insert) in the compact regime 
ranges from 1 to 2 mm (or approximately from 6 to 12 
bead diameters). In all cases, the particle fraction in the 
compact regime does not differ by more than 0.02 from 
the particle fraction Cmax in a static packing. 

Decreasing the relative humidity H of the air down 
to 40% or using "clean beads" with a smooth surface, 
changes significantly the behavior of the grain flow: in 
this continuous flow of grains is observed (its ve- 

locity is roughly constant with time and of the order of 
a few centimeters per second). In return, for "rougher 
beads" or relative humidity H > 50%, the flow is inter- 
mittent and may remain blocked for a large fraction of 
the time, particularly near the hopper (observations us- 
ing a SEM microscope display particles with a size of 
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FIG. 2: Spatio-temporal diagram of the grain flow: (a) at the 
top of the tube (100mm to 110 mm from the hopper), (b) at 
the bottom of the tube (850 mm to 860 mm from the hopper). 
Each diagram corresponds to a time lapse of 0.6 s and to a 
field of view of 10 mm. 



up to several micrometers at the surface of the "rough 
beads" with spacings increasing with the size from a few 
/im up to a few 10 /im). Experiments reported here have 
been performed at a relative humidity H = (55 ± 5)% for 
which the grain flow is intermittent. At the outlet of the 
constriction, flow is modulated but does not stop com- 
pletely. The modulation may be periodic, particularly at 
large flow rates. In the remainder of the paper, we shall 
focus on this latter intermittent compact flow regime. 



A. Spatio-temporal characteristics of the 
intermittent flow 

Figure |5] displays two spatio-temporal diagrams of the 
grain flow over two 10 mm high segments located in the 
top and bottom parts of the tube. Time corresponds to 
the horizontal scale. The two diagrams were obtained at 
slightly different times during a same experiment: due 
to the excellent periodicity of the flow, they are however 
comparable. Time intervals during which grains are at 
rest in the tube are marked by horizontal line segments 
with constant grey levels; those during which grains are 
moving down are marked by inclined striations: the slope 
of which indicates the grain velocity. Starting from rest, 
the grains flrst accelerate and reach a constant velocity 
before flnally stopping abruptly. Figures HJ^a) and (6) 
show that the fraction of the total time during which 
grains are flowing is much larger at the bottom of the 
tube than at the top. 

A quantitative comparison is made in Fig. in which 
the duration Tf of the flow phase at two different dis- 
tances from the hopper and the period T of the process 
are displayed as a function of the time average q of the 
superflcial velocity (T is the sum of the durations of the 
flow (Tf) and static phases). The time lapse Tf is 40% 
lower in the upper section than in the lower one; it in- 
creases linearly with q in both cases, while, on the con- 
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FIG. 3: Duration of flow phase T/, as a function of the time 
average q of the superficial velocity, at two different distances 
z from the top of the tube: z = 250mm (A) and 2 = 600 mm 
(a). Global period of the intermittent flow T (■). 
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FIG. 4: Mean grain velocity v vs. mean superflcial velocity 
q for 2 = 250 mm (T, left axis) and z = 600mm (V, left axis). 
Displacement Sz of the grains during a flow period as a func- 
tion of g (■, right axis). 
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FIG. 5: Spatio-temporal diagram of the grain flow at dis- 
tances between mm and 50 mm from the hopper at the top 
of the tube (g = 0.07 ms~^); total time lapse corresponding to 
diagram: 0.7 s. 

on Fig. 21 For large mean flow rates, the instantaneous 
grain flow rate q in the tube exceeds the maximum flow 
rate from the hopper. A low particle fraction bubble thus 
builds up at the top of the tube: it appears as a light 
zone in the spatio-temporal diagram of Fig. |S1 (the local 
particle fraction in this bubble may be lower than 0.2). 
The abrupt initial downwards slope of the contour of the 
light zone reflects the fast downwards motion of the upper 
boundary of the compact packing below the bubble. As 
soon as the motion of the packing stops, the bubble gets 
filled up from the hopper, although at a lower velocity 
that decreases with time (this may be due to the build 
up of an adverse pressure gradient between the bubble 
and the hopper). 




trary the period T decreases with q. Spatio-temporal 
diagrams also allow one to estimate qualitatively the ve- 
locity of the grains in the flow phase from the slope of the 
striations. Figures|2Ia) and (&) show that the grain speed 
is much higher at the top of the tube than at the bottom 
as required by mass conservation. The particle fraction 
c varies little along the tube and the global displacement 
of the particles during one period of the flow must be 
the same at all heights. The mean value v of this veloc- 
ity during the flow phase can be estimated quantitatively 
from the relation: 



The variation of v with q in both sections of interest 
is displayed in Fig. 01 in both cases, v increases roughly 
linearly with q. The displacement 6z of the grains during 
one flow period is equal to v x Tf and is also plotted 



B. Spatial variations of the intermittency effect 

Information obtained from the spatio-temporal dia- 
grams are complemented by measurements of the local 
pressure variations. Figure El displays time recordings 
of the pressure variations Sp = p — Po measured on four 
sensors at different heights z from the top {po is the atmo- 
spheric pressure). These curves are nearly periodic; the 
pressure drops sharply while grains are flowing (by up 
to 3000 Pa near the top of the tube). Then the pressure 
increases back above Po and decays slowly (or remains 
constant) while the grains are at rest. The time average 
of the pressure is close to zero (at most 25 Pa) for all four 
transducers used for the measurements. 

The flow has been filmed in the intermittent regime us- 
ing a 1000 frames per second high speed camera equipped 
with analog input channels connected to the pressure 
transducers. The onset of the fiow at a given height is 
observed to coincide exactly with the beginning of the 
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FIG. 6: Time recording of pressure variations 5p at four dif- 
ferent distances z from hopper (from top to bottom, z — 200, 
450, 700 and 950 mm), for g = 0.02 ms"\ Dashed lines corre- 
spond to the onset and the blockage of the flow. 
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FIG. 7: Processed spatio-temporal diagram over a 500 mm 
length of tube (200 mm < z < 700 mm) during one period of 
the intermittent flow. Light zones correspond to the static 
phase and the darker region to the flow. 



pressure drop; the transition is less sharp for the block- 
age since pressure diffusion is still significant even af- 
ter the grain flow has stopped. The region between the 
two dashed lines in Fig. |B1 corresponds to the flow phase 
(the lines are drawn, within experimental error, between 
points corresponding to times at which the flow is ob- 
served to start or stop). Comparing the widths of this 
region in the different curves confirms that the fraction 
of time corresponding to the fiow phase increases with 
z. On the contrary, the amplitude of the pressure drop 
decreases with z. 
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FIG. 8: Variation of the absolute velocities of propagation of 
the decompaction \vdf\ (▼) and blockage \vtf \ (■) fronts with 
the mean superficial grain velocity q. 



The propagation of the onset and blockage of the flow 
were studied more precisely on spatio-temporal diagrams 
of a 500 mm high region of the tube (Fig. d). The dia- 
gram has been processed by subtracting each linear im- 
age from the next one to detect more precisely the onset 
and blockage of the flow. Domains corresponding to the 
static phase have a light shade (the difference is not zero 
because of the noise of the camera and the fluctuations of 
the illumination). The domain associated with the flow 
phase is globally darker (the fluctuations of the grey lev- 
els are large because one subtracts specular reflections of 
light on the beads varying very much during the motion). 
The velocity \vdf\ for the onset of the flow corresponds to 
the slope of the left boundary in Fig. [7| it was roughly 
constant in that experiment but decreases slightly with 
distance in others. The velocity |uf,/| for the blockage is 
given by the slope of the right boundary (averaged along 
the tube) and is lower than \vdf\- Both the onset and 
the blockage propagate upwards. The variations of \vbf\ 
and \vdf\ with q are plotted in Fig. |H| The velocity \vbf\ 
for the blockage increases linearly with q and ranges be- 
tween 20 and 60 m s~^ . This variation will be discussed in 
section lTV Al On the contrary, the velocity \viif\ of the lo- 
cation of the onset of the grain flow (decompaction front) 
is independent of the mean superficial velocity and is of 
the order of 11 ms^^. This implies that the decompaction 
process is independent of the downwards structure of the 
fiow and of the width of the constriction. 

This velocity value may be related to the velocity of 
sound in two-phase systems with a large density and com- 
pressibility contrast between the two phases (water-air 
bubbles for instance). In this case, by assuming an isen- 
tropic expansion of the gas, one may compute the velocity 
of sound Wsound via 18]: 
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pc(l - c) ■ 



(4) 



in which c and p are the volume fraction and the den- 
sity of the dense phase, Po is the atmospheric pressure 
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and 7 = 1.4 for air which is considered as having neghgi- 
ble density. This amounts to considering the air-grain 
system as a homogeneous mixture of both high com- 
pressibiUty and density resulting in a low sound veloc- 
ity. Taking p to be equal to the density of the grains 
and c to the mean particle fraction c ~ 0.615 leads to 
I'sound — 15ms~^ which is indeed of the order of the 
value \vdf\ — 11 ms~^ deduced experimentally (0.615 rep- 
resents a typical value of the average particle fraction as 
estimated from the measurements displayed in Fig. [TUI 
assuming as above that Cmax = 0.63). 

Data points from the four curves in Fig. Altogether with 
the boundary conditions p = Po at both ends of the tube 
provide the overall shape of the pressure profile along the 
tube. The spatial variation of Sp = p — po at several dif- 
ferent times during one same period of variation of the 
flow within the periodic sequence of Fig. El is displayed in 
Fig. El The thick solid line in Fig.jSJa) corresponds to the 
static phase. A pressure maximum subsists in the cen- 
tral part of the tube from the previous flow cycle and the 
associated pressure gradients drive air towards both ends 
of the tube through the motionless packing. When flow 
is initiated, a low pressure zone appears at the bottom 
of the tube where grains are moving, and propagates up- 
wards. The amplitude of the pressure minimum increases 
strongly during this propagation. This amplification ef- 
fect is discussed below in section [3 After the pressure 
minimum has reached the top of the tube, pressure starts 
to increase again [Fig. (5^6)]. At first, the local particle 
fraction reaches a minimum and the expansion of the gas 
stops. Then, when the granular fiow stops in the tube 
(first at the bottom) , the particle fraction reverts back to 
its static value. Thus, air gets compressed above the at- 
mospheric pressure since some additional air has leaked 
between the grains during the flow period (see for in- 
stance the curve corresponding to i = 172 ms in Fig.jH^b)). 
Finally, the pressure profile along the static packing re- 
laxes back to its initial distribution at the beginning of 
the flow phase and a new fiow cycle may start again. 

Note that when the intermittent fiow regime does not 
appear spontaneously, it can also be triggered externally 
by periodic injections and suctions of air close to the 
constriction at a frequency comparable to that of the 
natural oscillations. The other characteristics of the fiow 
such as the amplification of pressure and particle fraction 
variations as they propagate upwards are the same as in 
the spontaneous intermittent flow regime. 



C. Particle fraction and flow rate variations 

Spatio-temporal diagrams have allowed us to identify 
granular motions in the various phases of the intermit- 
tency process; their relation with the spatial distribution 
of the pressure 5p has also been analyzed. We proceed by 
discussing fiow rate and particle fraction measurements 
which provide additional information on the fiow struc- 
ture. 
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FIG. 9: Pressure deviation {5p = p — po) variations with 
distance along the tube for several measurement times, for 
a mean superficial velocity g = 0.02ms~'^. (a) Sofid line 
obtained during the static phase at t = s. Dashed lines 
are profiles obtained at various times during the flowing 
phase: t=16ms, t = 31ms, t = Q2ms, t = 86ms, t= 109ms 
and t = 133 ms (dash length and spacing increase with time) . 
(fe) Solid line obtained near the end of the flowing phase 
at t = 133 ms. Dashed lines are profiles obtained at var- 
ious times: t — 141 ms, t = 148 ms, t — 156 ms, t = 172 ms, 
t — 211 ms, t = 266 ms (dash length and spacing increase with 
time) . 



The variations of the particle fraction c and pressure 
8p at a fixed height, and of the fiux Q^^'^^ of air into the 
hopper are displayed in Fig. ^1 Variations of the local 
particle fraction c associated with the onset and blockage 
of the fiow are only of the order of 0.015 in the present ex- 
ample and were never higher than 0.02. Moreover, these 
variations are evidently synchronous with those of the 
local pressure p in the fiowing phase. Note the slow pres- 
sure decay when the grains are at rest. The amplitudes of 
the variations of c and bp decrease, while their duration 
increases with distance from the hopper. In addition to 
demonstrate this close correlation between the variations 
of p and c, our study indicates an increase of the am- 
plitude of the c variations with the instantaneous grain 
velocities. 

The time recording of Q^'^'^ indicates a strong increase 
of the downwards air fiow rate during the fiow phase. 
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FIG. 10: Time recordings of the pressure deviations 5p and 
particle fraction c variations (at z = 200 mm) and of the mea- 
sured air flow-rate Q'^'"^" into the hopper. Q™"'''' is defined 
as positive for air flowing downwards. The dashed line corre- 
sponds to the onset of the flow. 



This increase reflects the entrainment of air by the gran- 
ular flow and the compensation of the volume of grains 
leaving the hopper by the inflow of air. The variation of 
Qmeas slightly delayed relative to those of p and c since 
a time is required for the grain flow to propagate to the 
top of the tube and into the hopper. A transient upward 
gas flow {Qa^'^^^ < 0) is observed after grains stop moving 
in the upper part of the tube, corresponding to air being 
expelled through the grain packing. 



IV. QUANTITATIVE ANALYSIS OF THE 
INTERMITTENT FLOW COMPONENTS 

A. Conservation equations for air and grains 

In the following, z = corresponds to the top of the 
tube and the coordinate z increases downwards. Down- 
flow is defined as being positive. v{z,t) and u{z,t) are, 
respectively, the grain and air velocities (averaged over 
the flow section) in the laboratory frame at a time t and 
a distance z from the hopper. They are related to the 
superflcial velocities q{z,t) and qa{z,t) and to the local 
particle fraction c(z, t) (also averaged over a flow section) 
by the relations: 



cv, 



9a = (1 - c) u. 



(5) 



(6) 



The superficial velocity q of the grains is related to the 
variation of c with time by the volume conservation equa- 
tion: 



dc 
di 
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A similar relation can be written for the mass flux of air 
Paqa (Pa being its density): 



(8) 



Assume that the variations of the gas pressure are isen- 
tropic so that: 



P ^ 
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constant. 



(9) 



{po is the atmospheric pressure, pa„ the air density at po 
and 7=1.4). Then Eq. (|SJ| becomes: 
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dt 
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(10) 



A first prediction provided by mass conservation condi- 
tions is the velocity Vb / at which the blockage front moves 
up the tube. We define Cmax to be the particle fraction 
in the static column and c+ and the particle fraction 
and grain velocity right above the front. Vbf is obtained 
by expressing the conservation of mass of the grains in 
a slice of thickness 6z while the front moves from the 
bottom to the top of the slice: 



(11) 



Since increases linearly with the superflcial veloc- 
ity q, equation l(TT|l implies therefore that \vbf\ should 
also increase linearly with q. This agrees with the trend 
displayed in Fig. |S1 Using the experimental values of 
, and Cmax provides an order of magnitude of vi,/ 
{\vbf\ ~ 20ms~^) which is in a good agreement with that 
obtained in Fig.|Sl 



B. Relative magnitude of air flow components 

Another issue regarding the flow of air through the 
system is the relative importance of the passive drag of 
air by the moving grains and of the permeation of air 
through the grain packing. The first flow component is 
equal to q{l—c) jc and the second is assumed to be related 
to the pressure gradient by Darcy's equation: 



Darcy _ K dp 



(12) 



where rja is the viscosity of air and K the permeability 
of the grain packing. The superficial velocity qa of air is 
then given by: 



K dp q 
qa = T- + - 1 - c). 

T]a oz C 



(13) 



Equations lO, (|10|l and IjlSI) can be combined to yield: 
dp 7p dv 



dp 
di 



dz (1 — c) dz 



d^p , , 



0.00 0.05 0.10 0.15 

7// 7- (ms"') 

FIG. 11: Average inflow of air as a function of q/c. 



in which the coefficient D satisfies: 



D 



"fpK 

(1 - C)Va ' 



(15) 



(the second order term [D/{'-fp)]{dp/dz)'^ has been ne- 
glected in Eq. (|14|l 'l. Considering locally the system as a 
packing of spherical beads of diameter d, the permeabil- 
ity K of the grains can be estimated from the Carman- 
Kozeny relation: 



K = 



(1 



180 c2 



(16) 



In the following, we assume that to leading order, p — Po 
and c = c, so that D and K are constants. 

The average of c^^^^^ over a flow period is determined 
by plotting the time averaged value of Q™"^ as a function 
of the ratio g/c in Fig. ^2 (this amounts to assuming 
that the mean grain velocity is equal to q/c which is 
approximately verified due to the small relative variations 
of c in these experiments). Volume conservation implies 
that Q^^^^^ /ttR^ and the sum g + g^T of the volume flow 
rates of air and grains are equal [see Eq. ©J. Using 
Eq. 1^21), this leads to: 



<la 



Darcy 
Qa 



(17) 



Experimentally, the variation of Q^°^ with q/c is lin- 
ear with a slope equal to 1±0.02 in Fig. ^2 These re- 
sults indicate that the time average of q^^"''^ is negligi- 
ble within experimental error although its instantaneous 
value is non zero. The relative motion of air and of the 
grains is therefore significant only during transient flow. 



Pressure variation mechanisms 



0.605 



0.610 



0.615 



FIG. 12: Pressure vs. particle fraction variations during 10 
cycles of the intermittent flow. Arrows indicate the direction 
of the variation on each part of the curve: the vertical phase 
corresponds to the slow pressure relaxation through the mo- 
tionless grain packing, the lower section corresponds to the 
onset of the flow and the upper one to the blockage of the 
flow. The dashed line has the theoretical slope 7Po/(l — c). 



pressure. We first assume that variations of density pa 
are only due to those of the local particle fraction c. This 
amounts to neglecting the permeation of air relative to 
the grains. Then, dqa/dz and [qa/{ip)] dp/dz are equal 
to zero in which case Eq. (|10|l becomes: 



1 



dc 



(1 



dt 



1 dp 
IPo dt ' 



(18) 



(j) and c are replaced by Po and c within a first order ap- 
proximation) . This relation is tested in Fig. El by plot- 
ting pressure variations as a function of particle fraction 
variations: for variations small compared to the mean 
value, data points should be located on the dashed line 
of slope 7Po/(l — c) shown on Fig. 1121 The two param- 
eters are clearly strongly correlated and the slope of the 
curves has the right order of magnitude in the initial and 
final phases of the variation. However, pressure is higher 
in the blockage phase than during the development of 
the fiow, indicating that dqa/dz cannot be neglected in 
Eq. ifTUI) (the last term [qa/[ip)\ dp/dz is of lower order 
and should remain negligible). Some air gets sucked into 
the lower pressure regions when grains flow; this leads to 
an overpressure when flow stops and the particle fraction 
increases. 

This overpressure decays slowly during the remainder 
of the static phase (short vertical part of the curve at top 
right) due to the finite permeability of the packing. In 
this static phase of the flow sequence, u = 0. The two 
first terms on the right hand side of Eq. 114|) are equal 
to zero leading to the diffusion equation: 



dp 
'dt 



D 



d'^p 



(19) 



The dynamical properties of the system are largely de- 
termined by spatial and temporal variations of the air 



Experimentally, the particle fraction varies from 0.605 
to 0.625; the permeability K computed from Eq. H16|) 
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ranges therefore from 2310-^2 2910-^2m-2. Tak- 
ing rja = 1.85 10~^ s~^, Po = 10^ Pa and a mean value 
of K (iC^ 26 10^12 one obtains D^O-Sm^s"!. 

The characteristic diffusion distance during the period 
T ~ 0.25 s of the intermittent flow is ^/DT ~ 0.35 m 
which is comparable to the length over which pressure 
fluctuations diffuse. During the onset of the flow which 
takes place over a shorter time (a few hundredths of a 
second) the diffusion of the pressure variations is suffi- 
ciently small for Eq. IjlHI) to be approximately valid. One 
thus concludes that the permeability of the grain packing 
is high enough to allow for relative motions between air 
and the grains during the flow phase; but too low to al- 
low for a global mean flow of air through the packing. In 
the previous computation, the initial particle fraction of 
the static packing is taken equal to Cmax = 0.63. Assume 
that wall effects reduce this value to 0.61 as discussed in 
section m so that the particle fraction varies from 0.585 
to 0.605 during a flow cycle. The mean corresponding 
value of K is 32 10^^^ voT^ so that D would increase to 
~ 0.6 m^ s~^ and remain of the same order of magnitude 

as for Craax = 0.63. 

Our experiments suggest that the intermittent flows 
result from the amplification of pressure and velocity 
fluctuations appearing spontaneously or induced in the 
constriction at the bottom of the tube. One expects the 
period of the intermittency to be determined by the char- 
acteristics of the flow over the full length of the tube and 
not just at the constriction. The period of the flow is of 
the order of the characteristic diffusion time of air over 
a distance of half of the tube length. This diffusion time 
represents the order of magnitude of the time necessary 
for pressure perturbations induced during the flow period 
to relax before flow is reinitiated. It may thus be closely 
related to the duration of the static phase of the periodic 
flow. 



V. MODELING OF THE INTERMITTENT 
FLOW 

In this section, the amplification of the decompaction 
wave as it propagates towards the top of the flow tube is 
modeled numerically. This approach yields insight into 
the dependence of the amplification on the permeability 
of the grain packing. 



A. Equations of motion and force distributions 

The equations of motion and force distribution in the 
granular packing will first be written. They will then be 
solved numerically together with the conservation equa- 
tions for air and grains. For simplicity, one models the 
case of a lower particle fraction zone propagating upwards 
in a vertical channel. 

The equation of motion reflecting momentum conser- 
vation for both air and grains contained in a vertical slice 



with tube radius i? can be written as: 

dv dp 9(7^2 2 

(the equation has been divided by the section S — -kR^). 
Acceleration terms for air have been neglected due to the 
large density difference between the air and grains. The 
friction forces on the air applied at the walls are similarly 
negligible due to their low area compared to that of the 
grains. 

Note that when Eq. (excluding the gravity, Uzz, 
and a^r terms) is combined with Eq. H14|l (neglecting the 
second order v dp/dz and diffusion terms), one obtains a 
sound wave propagation equation with a velocity given 
by Eq. Q. 

The left side of Eq. (|20|l corresponds to the accelera- 
tion term and the first term on the right to the weight 
of the grains, azr and Ozz represent, respectively, fric- 
tion forces on the grains at the lateral walls and verti- 
cal stress forces between grains in the tube section; they 
are assumed to be constant around the perimeter of the 
tube and in each section. cr,r is then estimated by ap- 
plying Janssen's model 0, 0, 0| assuming first that 
fzr — lJ-(Jrr foUowing Coulomb's relation, a^r is taken 
to be Kjazz, where the coefficient Kj depends on the 
structure of the packing and of the surface properties of 
the grains. Equation l|20() then becomes: 

dv dv dp dcTzz 

in which the characteristic length A is given by: 



In a static packing, all terms involving the grain veloc- 
ity vanish: the weight of the grains is balanced by the 
pressure gradient and the stress forces at the walls and 
in the tube section. When flow starts, the local parti- 
cle fraction decreases [Fig- lior b)] which, in turn, reduces 
friction forces at the walls, resulting in an acceleration of 
the flow. No exact relation is available on the variations 
of the stress forces. It will be assumed in the following 
that the resultant of the weight of the grains, the pressure 
gradient and the friction forces increases linearly with the 
deviation of the particle fraction c from its static value 
Cmax- Equation H21|l then becomes: 

dv dv / 1 dp\ 

-TTT ^ -V— + J[Cmax- C) [g — , (23) 

dt dz \ pcdz J 

where J is the parameter characterizing the magnitude 
of the variation. This assumption represents a first order 
approximation: it will be valid in the initial phase of the 
development of the perturbation as long as the derivative 
of the resultant force with respect to the particle fraction 
is non zero for c = Cmax- Other types of functions may be 
envisioned to replace J{cmax — c) : many will also repro- 
duce an amplification phenomenon provided they predict 
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a resultant force in the right direction with a monotonous 
dependence on {cmax — c) and if (as will be seen below) 
the prefactor equivalent to J is large enough. 



B. Numerical procedure 

Numerical simulations are realized on a simplified ID 
model in which a region of reduced particle fraction 
(called bubble) moves upwards inside a vertical tube filled 
with an otherwise grain packing. The bubble length is 
much smaller than that of the tube and it is assumed 
to be far enough from cither of the ends that boundary 
conditions do not influence the motion. As suggested 
by the experimental observations, the particle fraction 
and the velocity are discontinuous at the bottom end of 
the bubble (the blockage front), while they vary continu- 
ously at the top (the decompaction region). Air pressure 
is continuous in both cases but the pressure gradient is 
discontinuous at the lower front. 

In order to solve the problem numerically, one uses the 
conservation equations for grains Q and air H10() estab- 
lished in section IIVI and the equation of motion H23|) of 
the grains. For the sake of practicality, Eq. (|14l) . in which 
Darcy's law is used to determine qa, replaces Eq. (|ll)l) . In 
the static regions of the packing, all terms in Eqs. ((TJ and 
(|23|l are equal to zero. Equation H14(l reduces then to the 
diffusion equation (|19|l for air with a coefhcient D given 
by Eq. (|T^. 

At large distances below and above the bubble, the 
air pressure p is equal to the atmospheric pressure po] 
c and V are respectively equal to Cmax and zero both 
far above the bubble and everywhere below the blockage 
front. Additional boundary conditions corresponding to 
the mass conservation for air and for the grains must 
also be verified at the lower front (in addition to pressure 
continuity). In the reference frame moving with the front 
at Vbf, grain mass conservation is expressed by Eq. 111|) 
which can be rewritten as: 

(u+ - t;6/)c+ = -VbfCrnax- (24) 

In the same way, the mass conservation of air implies; 



(v+-Vbf)il-c+)- 



K dp 

Va dz 



= -Vbf{l-Crnax)- 



K dp 
rja dz 



(25) 

The indexes + and correspond to values close to the 
blockage front respectively inside the bubble and in the 
static zone. Combining Eqs. (|24|l and (|25fl leads to: 



K dp 
rja dz 



K dp 
rja dz 



(26) 



The initial particle fraction, velocity and pressure pro- 
files are assumed to follow a l/cosh(2;) variation with a 
cutoff at the lower front (z = Zf) for c and v. Pressure 
is continuous at z = Zf with a faster decrease below. 



The corresponding profiles are displayed in Figs. llSf a- 
c) (rightmost solid curves). The initial particle fraction 
and pressure variations are related by Eq. Q above the 
front (this amounts to assuming that the bubble has 
been created through a local decompaction of the grains 
inducing an isentropic adiabatic expansion of the air). 
The amplitudes of these initial variations are respectively 
(5c = 0.001, (5p = 600Pa and (5u = 3 10"^ ^^g-i 

Pressure, particle fraction and grain velocity profiles at 
subsequent time steps are then computed by solving nu- 
merically Eqs. (|7|), (|14|l and H23() in the reference frame of 
the front using a finite difference Euler explicit scheme. 
For the resolution of Eq. (|23|) . an upwind scheme was 
required to avoid numerical instabilities due to the ad- 
vection term vdv/dz. The global length of the simula- 
tion domain is 1.4 m and the mesh size 2 mm. The origin 
of distances is taken at the initial location Zf{Q) of the 
front (the vertical dashed line on Figs. Wit a-c)). Simu- 
lations correspond to a time lapse of 27 ms during which 
the bubble propagates over a distance of 0.4 m. At each 
step of the simulation, the pressure, the particle fraction 
and the velocity are determined at all distances from the 
front and the front velocity is also obtained. This al- 
lows one in particular to determine the new location of 
the front after each time step and to translate distances 
from the front into distances in fixed coordinates which 
will be used in subsequent plots. 



C. Numerical results 

1. Qualitative observations 

Particle fraction, velocity and pressure profiles at dif- 
ferent times obtained for a typical set of values of J, K 
and Cmax are displayed in Fig. E| The other parameters 
of the simulation are given by the experimental values 
listed above or computed from J and K . A clear ampli- 
fication effect is observed for all three parameters. The 
velocity v reaches its maximum value faster than both 
the particle fraction c and the pressure 5p (typically af- 
ter 0.3 m versus 0.5 m). The ratio of the maximum and 
initial amplitudes is also higher (20 for v versus 2 and 
1.6 for 5p and c). A small localized overpressure appears 
in the static region close to the initial location of the 
bubble: it is due to air advected by the grains into this 
region at short times [from Eq. (|26|l ]. At later times, the 
pressure gradient close to the front increases so that the 
diffusion of air back into the bubble balances the advec- 
tion term. The increase of the amplitude of the variations 
of 5p and c is consistent with the decrease of the width 
of these curves: the total area under these curves must 
remain constant as required by mass conservation of air 
and grains. The velocity v^f of the lower front is also 
roughly constant with a value \vbf \ ~12.5ms~^ of the 
same order of magnitude as that observed in our experi- 
ments. 

While the computed amplitude of the pressure vari- 
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FIG. 13: Numerical pressure (a), particle fraction (b) and ve- 
locity (c) profiles, obtained for J = 450, /sT = 20 10"^^ and 
Cmaa; = 0.63. Froui right to left: solid lines: initial profile at 
t = s - dashed lines: t = 9 ms, t = 15 ms, t = 21 ms, t — 27 ms 
(dash length and spacing increases with time). The vertical 
dashed lines represents the position of the front Zf at t = 0s. 



ations is of comparable magnitude to the experimental 
ones, particle fraction and velocity variations are smaller. 
This is likely a result of the difference between the model 
case of a localized bubble and the actual intermittent 
flow which sometimes extends over the entire length of 
the tube. 



14- 


^^^^^^^^^^^^^^^^ 


-3000 


12- 




zuuu tr 


" 10- 
B 8- 


^^^^^^^^ 


*-< 
K 

- [000 — . 
IP 


6 - 


1 / 


tn 


> 4- 


Y 


-0 ^ 


2- 






n- 




--1000 



100 



200 



300 



400 



500 



FIG. 14: Variation of the amplification parameter Aylp/Az 
(o) and of the absolute front velocity (■) with the force 
parameter J, for K = 2Q 10"^^ m^. 



2. Quantitative dependence on the modelization parameters 

The key adjustable parameter in the model is the co- 
efficient J in Eq. H23|l (that characterizes the effective 
acceleration of the grains when the particle fraction c 
decreases below Cmax)- Figure displays the depen- 
dance on J of the front velocity and of the param- 
eter AAp/Az characterizing the amplification. AAp/Az 
is determined by plotting the minimum value Ap of the 
pressure variation curves [as shown in Fig. Il^f g')] as a 
function of the coordinate 2/ of the front for each curve. 
After a short transient, Ap increases quickly with Zf and 
the variation levels off thereafter. The amplification pa- 
rameter AAp/Az is taken equal to the maximum ab- 
solute value of the slope \dAp/dzf\ of this curve. The 
permeability K of the grain packing has been taken 
equal to the value K = 2010~^^ m^). Using the values 
of {cmax — c) deduced from the simulations, the maxi- 
mum value J = 500 in Fig. E| corresponds in Eq. (|23|l to 
an effective acceleration of the grains of the order of 5/3 
or lower. 

For J — 0, AAp/Az is negative and initial variations 
of 6p and c damp out as the bubbles propagate upwards. 
AAp/Az becomes equal to for J = 40 and increases 
steadily at higher values. The front velocity \vbf\ also 
increases with J in the amplification domain but never 
becomes zero. Its value is of the same order of magnitude 
as the experimental front velocities, particularly when J 
is large. We note that in the case of a bubble of slowly 
varying length, Vdf and Vbf are nearly equal while they 
differ in the experiments so that more detailed quanti- 
tative comparisons would be irrelevant. These results 
indicate that the simple model discussed above allows 
us to reproduce the amplification effect observed exper- 
imentally. In the following, it is applied to study the 
influence of the permeability K of the grain packing on 
the amplification process. 

The permeability influences the dynamics of the sys- 
tem through the diffusion coefficient D which is pro- 
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FIG. 15: Variation of the amplification parameter AAp/ /S.z 
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portional to K and appears in Eq. (|14(l . The depen- 
dance of \vbf \ and the amphfication parameter A^p/Az 
on K are displayed in Fig. 1151 The range of perme- 
abilities investigated is between 0.5 and 4 times the 
value if = 20 10^^^ m^ assumed in Fig. The front 
velocity viyf is first almost independent of K while 
^Ayl t^z sharply decreases with K and vanishes for 
K ~8010~^^m^. Note that reducing Cmax from 0.63 to 
0.61 assuming the influence of wall effects increases the 
permeability K from 22 to 27 10~^^ m^ and only reduces 
the amplification parameter by roughly 15%. 

The above results demonstrate that the diffusion term 
proportional to D (and therefore to K) in Eq. H14() is 
purely dissipative and acts only to enhance the attenu- 
ation of the decompaction wave. This results from the 
fact that, in high permeability packing, pressure varia- 
tions induced by the decompaction quickly diffuse away 
from the front. On the contrary, in low permeability sys- 
tems, these pressure variations remain localized: their 
gradients accelerate nearby grains and so sustain and am- 
plify the motion. Intermittency effects can therefore be 
expected to be stronger for flowing grains of smaller di- 
ameter. 



VI. CONCLUSION 

An experimental study has allowed us to analyze the 
dynamical properties of periodic intermittent vertical 
compact flows in a tube. Our experiments involved si- 
multaneous measurements of local pressure and particle 
fraction fields, combined with spatio-temporal diagrams. 

A key observation is the amplification of the parti- 
cle fraction and pressure variations as the decompaction 
front marking the onset of the flow propagates upwards. 
This propagation takes place at a much higher veloc- 
ity >10ms~^) than that of the grains (a few 
10~^ms~^). The velocity \vhf \ of the blockage front is 



still faster than \viif\: it increases with the flow rate and 
ranges between 20 and 50ms~^. 

The characteristic grain velocity during the flow phase 
increases with height above the outlet: this variation is 
associated with a decrease with height of the duration Tj 
of the flow phase as required by the mass conservation 
of the grains. Close to the top of the flow tube, Tf rep- 
resents only 10% of the total flow period T. The grain 
velocities and amplitude of the pressure fluctuations are 
sufficiently high in this zone that flow from the feeding 
hopper gets blocked and a transient low particle fraction 
region appears briefly at the top of the tube. 

A characteristic feature of these flows is the small am- 
plitude of the particle fraction variations (0.02-0.03 at 
most); this is much smaller than variations observed in 
the decompaction of a granular column with an open 
bottom end or a flow in the density wave regime. It is 
noteworthy that these small particle fraction variations 
nevertheless induce large pressure drops (up to 3000 Pa 
at the top of the tube). The decompaction process is 
indeed fast enough for air to undergo an isentropic ex- 
pansion: pressure variations closely mirror the variations 
of the local particle fraction. Diffusion of air through the 
grain packing is too slow to influence strongly that phase; 
however, it plays an important role in the longer static 
phase by allowing the decay of the pressure towards equi- 
librium before flow is reinitiated. The characteristic time 
for air diffusion through the packing is of the order of the 
total period of the intermittent flow and may be a key 
factor determining its value. 

The amplification effect can be reproduced by a ID 
model assuming an effective driving force on the grains 
that increases with the deviation of the particle fraction 
from the static value. This numerical model also indi- 
cates that the permeability K of the packing is a very 
important factor. The amplification effect (and perhaps 
also the intermittent ffows) should be observable only if 
K is sufficiently low that pressure perturbations do not 
diffuse away too quickly. 

Another important parameter is the relative humid- 
ity of the surrounding air which influence strongly the 
interaction forces both between the grains and between 
the grains and the tube walls. Preliminary results in- 
dicate that the intermittent flow regime is rarely ob- 
served at relative humidities of the order of 40% or be- 
low, at least with the beads used in the present study. 
Systematic studies in a broad range of relative humidi- 
ties should greatly help understand better the interaction 
forces which determine the characteristics of these inter- 
mittent flows. 
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